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Quantum mechanical superexchange interactions form the basis of quantum 
magnetism in strongly correlated electronic media. We report on the direct 
measurement of superexchange interactions with ultracold atoms in optical 
lattices. After preparing a spin-mixture of ultracold atoms in an antiferro- 
magnetically ordered state, we measure a coherent superexchange-mediated 
spin dynamics with coupling energies from 5 Hz up to 1 kHz. By dynamically 
modifying the potential bias between neighboring lattice sites, the magnitude 
and sign of the superexchange interaction can be controlled, thus allowing the 
system to be switched between antiferromagnetic or ferromagnetic spin inter- 
actions. We compare our findings to predictions of a two-site Bose-Hubbard 
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model and find very good agreement, but are also able to identify corrections 
which can be explained by the inclusion of direct nearest-neighbor interac- 
tions. 

Quantum spin systems on a lattice have served for decades as paradigms for condensed 
matter and statistical physics, elucidating fundamental properties of phase transitions and acting 
as models for the emergence of quantum magnetism in strongly correlated electronic media. 
In all these cases, the underlying systems rely on a spin-spin interaction between particles on 
neighboring lattice sites, such as in the Ising or Heisenberg model (|7]|2]|J]). As initially proposed 
for electrons by Dirac (0|5|) and Heisenberg (|2][6]), effective spin-spin interactions can arise 
due to the interplay between the spin-independent Coulomb repulsion and exchange symmetry 
and do not require any direct coupling between the spins of the particles. The nature of such 
spin-exchange interactions is typically short-ranged, as it is governed by the wave function 
overlap of the underlying electronic orbitals. In several topical insulators, such as ionic solids 
like e.g. CuO and MnO, however, antiferromagnetic order arises even though the wave function 
overlap between the magnetic ions is practically zero. In this case a "superexchange" interaction 
mediated by higher order virtual hopping processes can be effective over large distance 
which leads to an (anti) -ferromagnetic coupling between bosons (fermions) on neighboring 
lattice sites Q. Such superexchange interactions are believed to play an important role in 
the context of high-Tc superconductivity (|9l). Furthermore, they can form the basis for the 
generation of robust quantum gates similar to recent work in electronic double quantum dot 



systems (10 11 ), and can be used for the efficient generation of multi-particle entangled states 



(12 13), as well as for the production of many-body quantum phases with topological order 



q^\T5\\16^, 

We report on the direct observation of superexchange interactions with ultracold atoms in 



optical lattices (17 18). Previous experiments have shown that spin-spin interactions between 



neighboring atoms can be implemented in discrete time steps (|79] |20|) by bringing the atoms 
together on a single site and carrying out controlled collisions (|27]|2^|22]) or onsite exchange 



interactions ( [23| ). The superexchange interactions demonstrated here, however, directly imple- 
ment nearest-neighbor spin interactions in the many-body system and allow for a continuous 
"analog" simulation of spin lattice Hamiltonians. 



We probe the superexchange interactions by first preparing two atomic spin states of ^^Rb 
in an antiferromagnetic order (|2?1) and then recording the time evolution of the spins of neigh- 
boring atoms in isolated double well potentials (|25||2^|2^ for weak to strong tunnel couplings. 



For dominating onsite interactions over the tunnel coupling between lattice sites, we find pro- 
nounced sinusoidal spin oscillations due to an effective Heisenberg-type superexchange Hamil- 
tonian, whereas for stronger tunnel coupling a more complex dynamics emerges. In addition, 
we show how the strength and sign of the superexchange interaction can be directly controlled 
by introducing a potential bias between neighboring wells. Furthermore, we find that correc- 
tions to the two-site Bose-Hubbard model (BHM) which take into account the direct interaction 
between particles on neighboring lattice sites are needed to fully explain our data. 

Theoretical model. An isolated system of two coupled potential wells constitutes the simplest 
concept for the investigation of superexchange-mediated spin dynamics between neighboring 
atoms. We consider a single double well potential occupied by a pair of bosonic atoms with 
two different spin-states denoted by |t) and ||). If the vibrational level splitting in each well 
is much larger than all other relevant energy scales and intersite interactions are neglected, 
the system can be described in a two-mode approximation by a two-site version of the Bose- 
Hubbard Hamiltonian 



+ U (ntL^iL + ^tr^^r) , (1) 



where the operators and acrL,R create and annihilate an atom with spin a in the left and 
right well respectively, ncrL,R count the number of atoms per spin state and well, J is the tunnel 
matrix element, A the potential bias or tilt along the double well axis and U = U^y = g x 
J wl^^{-x)cPx the onsite interaction energy between two atoms in | |) and | |). Here, g = 
(47r/i^a|^)/M is the effective interaction strength with aj^ being the positive scattering length 
for the spin states used in the experiment, M is the mass of a single atom and t(;L,R(x) denote the 
wave functions for a particle localized on the left or right side of the double well. The state of 
the system can be described as a superposition of the Fock states {|t, |), ||, j), |t|, 0), |0, 11)}^ 
where the left and right side in the notation represent the occupation of the left and right well, 
respectively, and the states | jj, 0) and |0, t|) are spin triplet states. In the following, we will 



focus on the dynamical evolution of the population imbalance x = rii^ — and the Neel 
order parameter or "spin imbalance" A^^ = (n^L + ^^r — ^tR ~ ^iL)/2 starting with double 
wells initially prepared in |T, i). Here n^^j;L,R = (^ta;L,R) denote the corresponding quantum 
mechanical expectation values and nL,R = ti^l.r + ^il^r- 

In the limit of dominating interactions (U J), when starting in the subspace of singly 
occupied wells spanned by | T, i) and | j, T), the energetically high lying states | Ti, 0) and 
|0, tl) can only be reached as "virtual" intermediate states in second order tunneling processes. 
Such processes lead to a non-local (super) spin-exchange interaction, which couples the states 
1 1, 1) and I i, T) (see Fig. lA). More generally, for an arbitrary spin configuration with equal 



interaction energies U^^ = = 1/^ (see ref. {28)), the second order hopping events are 
described by an isotropic Heisenberg-type effective spin Hamiltonian in the limit [/ :» J (|3][75 



i^eff = —"^Je^Si, • aSr = — Jex (^'l ^^j^ + ^'l — 2Jq^SIS^ , (2) 

where S^^^ = |T)(i|L,R, S^^r = U)(T|l,r and Sl^^ = (n^L,R - ^iL,R)/2 denote the corre- 
sponding spin operators of the system, with S'l,r = Sx ± iSy. The effective coupling strength 
Jex represents the superexchange and can readily be evaluated by perturbation theory up to 
quadratic order in the tunneling operator which yields Jgx = 2 

When a potential bias A > is applied, the degeneracy of the two intermediate states in the 
superexchange process is lifted (see Fig. 1 A). For J, A <C t/ this leads to a modification of the 
effective superexchange coupling with now Jex = J'^ /{U+A)+J'^ /{U-A) = 2J^U /{U^-A^) 
( [75| ). By tuning the bias to A > [/, it is possible to change the sign of Jex and therefore to switch 
between ferromagnetic and antiferromagnetic superexchange interactions. For J <C |t/ — A| 
the picture of an effective coupling via two virtual intermediate states is again valid and the full 
reversal to Jex = — 2 J^/t/ is found to be reached for A = \/2 U . 

For symmetric double wells (A = 0), the Hamiltonian Eq. 1 can be diagonalized analytically 
to give a valid picture for all values of J and U within the single band BHM. A convenient 
basis is given by the spin triplet and singlet state \t/s) = (| |, |) ± 1 1, |))/v^ and the states 
|±) = (III, 0) ± |0, ||))/\/2. Two of the eigenstates are linear combinations of \t) and |+), 
where the one having the larger overlap with \t) is the ground state. The spin singlet 1^) and 
the state |— ) are already eigenstates themselves with energy and U respectively (see Fig. IB). 
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As a direct consequence, |— ) cannot be reached from the initial state |t, |) = (l^) + \t))/\/2. 
Therefore, the dynamical evolution of the spin imbalance contains only two frequencies 




^^i,2 = ^h/(^)' + i±i| . (3) 

The extraction of these frequencies from time-resolved measurements allows for the determi- 
nation of 2 J = h-yjuj^^ and U = h{ijJi — 002) within the BHM. 

As these frequencies can be measured with high accuracy, we are able to observe devia- 
tions from the simple BHM. We obtain a first correction by the inclusion of nearest-neighbor 
interactions (iJl) in an extended two-site Bose-Hubbard model (EBHM, see Eq. SI in ref. ([J/])). 
This modification introduces the inter- well interaction energy [/lr = 9^1 w'1{-k)w^{-x)(Px 
and a correction to the tunneling matrix element, which becomes = J + A J, where AJ = 
—g X J wI{-k)w^{x.)(Px. The inter- well interaction leads to a direct spin-exchange term, which 
in the limit U ^ J reduces the corrected superexchange coupling to J^^ = 2J'^ /U — [/lr- 
While we find that the corrections to the pure two-site Bose-Hubbard model are not negligible 
in the experimentally relevant parameter region, numerical calculations based on the multiband 
Schrodinger equation show that the direct exchange can never overcome the superexchange 
coupling term (see Fig. SI) and therefore change the nature of the ground state to be anti- 



ferromagnetic. This is in agreement with the Lieb-Mattis theorem ( p2| ), which states that the 
groundstate for two bosons has to be a spin-triplet state. 

Initial state preparation. In order to investigate the spin dynamics between neighboring 
atoms, we initially prepare a sample of ultracold neutral atoms with two relevant internal states 
It) and II) in a 3D array of double wells with Neel-type antiferromagnetic order iTiTiTi • • •) 
along one spatial direction (see Fig. 2A). State preparation was started by loading a ^^Rb Bose- 
Einstein condensate of typically around 8 x 10^ atoms in the |F = I.tyif = —1) Zeeman 
sublevel with no discernible thermal fraction from a magnetic trap with high offset field into 
a 3D optical lattice of double well potentials ( |25| ). This "superlattice" potential is obtained by 
superimposing on one axis two standing light fields with periodicity 382.5 nm (short lattice) 
and 765 nm (long lattice) and additional standing waves with periodicity 420 nm on the two 



perpendicular axes ( [27| ). Controlling all depths and the relative phase of the short and long 
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lattice allows to tune the double well configuration in terms of the Hamiltonian parameters 
J, U and A. The depths of the lattices are given in units of the short-lattice recoil energy 
Er = / (2MA^) with A = 765 nm throughout the article. The loading ramps were optimized 
to favor an occupation of two atoms per double well and to avoid heating to higher vibrational 
levels (31). After merging the double wells by ramping down the short lattice, a microwave 
rapid adiabatic passage was used to transfer all atoms into the |F = 1, = 0) state. Sub- 
sequently the magnetic trap is switched off and while keeping a homogeneous offset field of 
around 1.2 G, atom pairs were coherently transferred from \mF = Q'.rriF = Q) into spin triplet 
pairs (1+1;— 1) + | — 1; +1))/V^ by means of spin-changing collisions {33 34), The two mag- 
netic sublevels \mF = ±1) correspond to the two spin states |t) and ||). The remaining atoms 
in the \mF = 0) state, e.g. on singly occupied sites, are transferred into \F = 2,mF = and 
removed in a filtering sequence before the detection ^3l\. 

Finally, the short lattice was ramped up slowly in 20 ms thereby inhibiting a coherent split- 
ting of the atom pairs and leaving the double wells in a state with one atom on each side (2^. 
For the time of the ramp up, a magnetic field gradient of ^ 17G/cm in the direction of 
the superlattice is switched on. Therefore, the degeneracy of the states |t, |) and ||, j) in the 
double well is lifted by approximately 900 Hz, which enables an adiabatic loading of the state 
11,1) during the splitting process (see Fig. 2B). Numerical integration of a multiband ansatz 
for this procedure yields an expected fidelity of > 99% for creating an antiferromagnetic order 
along the axis of the superlattice. The mean population imbalance x{t) and spin imbalance 



Nz{t) of the ensemble of double wells was detected by applying a mapping technique (26 27) 
combined with a Stern-Gerlach filter (see Fig. 2C). A maximum spin imbalance of 60-70% was 
observed for our initial state corresponding to a probability of 80-85% for having prepared the 
desired state | T, i). We believe that this measured value is mainly reduced due to our detec- 
tion method. Direct spin-exchange processes emerge during the mapping sequence (23) and 
can lead to a mixing of the spin configuration and thus a reduction of the measured Neel order 
parameter (31). 

Time-resolved observation of superexchange interactions. The spin dynamics are initiated 
by rapidly ramping down the short lattice and thereby the double well barrier in 200 /xs, thus 
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significantly increasing the tunneling and superexchange couplings. After letting the system 
evolve for a hold time t, the spin-configuration was frozen out by ramping up the barrier in 
200 /xs, quenching both J and Jgx again. The measurement of the ensemble averages x{t) and 
Nz{t) is carried out as described above. 

Three typical time traces obtained by this procedure are shown in Figure 3. For low barrier 
depths (J/U > 1), we observe a pronounced time evolution of the spin imbalance Nz{t) con- 
sisting of two frequency components with comparable amplitudes and frequencies (Fig. 3A). 
With increasing interaction energy U relative to J, the frequency ratio increases, leaving a slow 
component with almost full amplitude and an additional high-frequency modulation with small 
amplitude (Fig. 3B). The fast component corresponds to first order tunneling due to the coupling 
of \t) and |+), which becomes more and more off-resonant as the barrier height is increased and 
therefore J/U is decreased. For J/U <C 1, it is completely suppressed and the only process vis- 
ible is the superexchange oscillation (Fig. 3C). For all barrier heights, the population imbalance 
x{t) stays flat, emphasizing that even though strong spin currents are present in the system, no 
net mass flow can be observed for our initial state. We fit the traces for Nz{t) with a sum of two 
damped sine waves with variable frequencies 001^2/ i"^^) and amplitudes 2. For the damping 
we assume a Gaussian characteristics with 1 / e-damping constants 71^2- The results of the fit are 
displayed in Fig. 4. For 1/short > 15 (inset in Fig. 4 A) we can identify only a single frequency 
component corresponding to the superexchange oscillation with 4 J^/ (hU) and full relative am- 
plitude (see Fig. 4B). We are able to observe this frequency down to 4.8(4) Hz at J/U = 0.023 
for l/short = 20 Er. The damping of the signal can be explained by the inhomogeneous distri- 
bution of coupling parameters due to the Gaussian shape of the lattice beams, which leads to a 
dephasing of the evolution within the ensemble. For l^hort ^ 17 E^, additional damping mech- 
anisms like tunneling to empty adjacent lattice sites (defects) or small residual inhomogeneous 
magnetic field gradients become relevant and limit the measurements (see Fig. 4C). 

The comparison of the results with the theoretical predictions by the simple BHM shows 
significant deviations at low barrier heights which cannot be explained by our uncertainties 
in the lattice depths. In this region, the EBHM can model the experimental data much more 
accurately. This can be understood by the fact that the inter- well interaction energy increases 
with decreasing barrier and begins to noticeably influence the dynamics i ^3l) . In fact, the EBHM 
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description yields h{uji — UJ2) ^ U + S^/lr and therefore directly explains the upward trend of 
this frequency difference for small short-lattice depths (see Fig. S2). For large barrier heights 
our experimental data is compatible with both models within the uncertainties of the lattice 
depths. However, here, the predictions of the EBHM are always closer to the measured values. 

Sign reversal of the effective coupling parameter. In order to demonstrate the controllability 
of superexchange interactions, the spin dynamics was investigated with an applied bias on the 
double wells for a short lattice depth of 15 and the same depths for the long and transverse 
lattice as before. Starting with an initial antiferromagnetically ordered state, as above, we first 
let the system evolve in symmetric double wells (A = 0) for to = 4.5 ms until the first node 
Nz{t) = of the spin imbalance is reached for the state {\t) + i\s))/^/2 (see Fig. 5). After 
freezing out the relative phase between 1^) and \t) by ramping up the potential barrier, a defined 
potential bias A is applied and a second evolution sequence with hold time f ^ t — to is 
initialized by ramping down the short lattice again to 15 Er. The subsequent detection follows 
the scheme described above. 

Fig. 5 A shows the evolution of the spin imbalance Nz{t) in symmetric double wells together 
with the time traces for two different bias energies A > U yielding an effective coupling of 
Jex(A) ^ — Jex(A = 0) and — Jex(A = 0)/2 respectively. The sign reversal of Jgx due to the 
introduction of the bias is directly visible by the change in slope of the spin imbalance t ^ to. 
It should be noted that the now negative sign of Jgx for bosons does not imply a violation of the 
Lieb-Mattis theorem, as the new ground state in this regime is the spin triplet state |ti, 0) and 
the superexchange couples the first and second excited states 1^) and \t) which have reversed 
order for A > f/ (see left inset in Fig. 5B). 

The introduction of a non-zero tilt leads to an increased sensitivity of the exchange fre- 
quency to fluctuations due to the inhomogeneities in the array of double wells which are most 
effective around A ^ U. Therefore the damping of the signal due to dephasing is stronger as 
A approaches U. A fit of a single dampened sine wave to the time traces obtained for various 
tilts yields the frequency curve shown in Fig. 5B together with the amplitude of the oscillation. 
Starting around 50 Hz, the oscillation frequency reaches a resonance for A/C/ ^ 1, where the 
amplitude reverses sign, leading to the observed time reversal in the dynamics. 
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Summary and outlook. We have demonstrated time-resolved measurements of superexchange 
spin interactions between ultracold atoms on neighboring lattice sites and have shown how to 
control such interactions with optical superlattices. Comparing the measurements to theoret- 
ical predictions of these spin interactions from first principles we find excellent agreement of 
our data to an extended two-site version of the Bose-Hubbard model. Although superexchange 
interactions become exponentially suppressed for deep optical lattices, the coupling strength 
2Jex/^ can be several hundred Hertz for lattice depths of around 12 — 15 and thus almost 
a factor of thousand larger than the direct magnetic dipole-dipole interaction of Rb atoms on 
neighboring lattice sites. One order of magnitude larger coupling strengths than the ones shown 
here, however, could still be achieved using electric dipole-dipole mediated spin interactions 



between ground state polar molecules ( p5| ). 

The demonstrated scheme to change the superexchange coupling strength and reverse the 
sign of the spin interaction can also be applied to the full ID chain, offering novel possibilities 
for engineering spin-spin interactions in optical lattices. It is now e.g. conceivable to engi- 
neer a setup with ferromagnetic interactions along one and antiferromagnetic interactions along 
another lattice direction. Furthermore, one can dynamically switch between ferro- or antifer- 
romagnetic interactions along a given lattice direction and follow the subsequent dynamical 
evolution of the quantum spin system. 

When the presented loading scheme is carried out without any magnetic gradient field dur- 
ing the splitting process, a valence-bond solid (VBS) type spin state can be efficiently 



engineered. Such VBS states can be viewed as a large array of robust Bell pairs ([J^lJS]). In 
principle, the superexchange interaction can be changed to be of Ising-type, e.g. by tuning 
the interspecies scattering length (T5\. Thereby, it can be used to create large entangled states 
out of the initially disconnected pairs, which have been shown to be powerful resources for 



measurement based quantum computation ([75]|59]). Moreover, controlling the superexchange 
interactions along different lattice directions also offers novel possibilities for the generation of 
topological many-body states for quantum information processing ( [74]|75| ). 
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Figure 1: Schematics of superexchange interactions. A Second order hopping via |Ti, 0) and 
|0, tl) mediates the spin-spin interactions between atoms on different sides of the double well. 
B Energy levels for A = and U ^ J. The evolution in the upper doublet of states (i) 
corresponds to the correlated tunneling of atom pairs (27), while the superexchange takes place 
in the lower one (ii). Both doublets are coupled by first order tunneling processes (iii). 
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Figure 2: State preparation and detection. A Spin triplet pairs are created on doubly occupied 
lattice sites and subsequently split under the influence of a magnetic field gradient to obtain 
antiferromagnetic Neel order. B Evolution of the eigenenergies with respect to the ground 
state energy during the splitting with V\ong = 10 E"^, Krans = 25 E"^ and a gradient of B' = 
17G/cm. The notation [^i, ^2] denotes the number of vibrational excitations for the first and 
second particle. C Detection of the population and spin imbalance. The population of the 
left well is transferred to a higher vibrational level of the underlying long-lattice well ( |25|[27| ). 
Subsequent band mapping and a Stern-Gerlach filter allow to determine x{t) and N^it) from 
the time-of-flight (TOF) images. 
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Figure 3: Spin and population dynamics in symmetric double wells. The time evolution of 
the mean spin Nz{t) (blue circles) and population imbalance x{t) (brown circles) are shown 
for three barrier depths within the double well potential: (A) ^short = QE^.J/U = 1.25, (B) 
^short = 11 Er, J/U = 0.26 and (C) l/short = 17 ^r, J/U = 0.048. The measured traces for the 
spin imbalance are fitted with the sum of two damped sine waves (blue lines). The population 
imbalance x{t) stays flat for all traces. 
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^short i^r) ^short i^r) 

Figure 4: Frequency components of the spin oscillations. A The frequencies obtained by fitting 
the spin imbalance data for various values of l/short (red and blue circles) as well as the frequency 
difference and the geometrical mean (green and black circles) are compared to the theoretical 
predictions of the simple BHM (colored regions). The width of the regions represents a 2% 
uncertainty in all lattice depths. The 2a error bars on the fit results and extracted values are 
in all cases smaller than the data points. The dashed lines are the predictions of the EBHM. B 
The plot of relative amplitudes vs. short lattice depth shows the suppression of the first order 
tunneling process towards small values of J/U. The damping coefficients for the two oscillatory 
components vs. lattice depth are displayed in C. 
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Figure 5: Sign reversal of the superexchange coupling. A After a quarter of a superexchange 
oscillation in symmetric double wells (black and gray dots, J/U = 0.08), a defined bias is 
applied, leading to a reversal in the slope of the spin imbalance and ongoing oscillations. For 
A/[/ = 1.32 (blue circles), the magnitude of the coupling is almost the same as for A = 0, 
whereas it is decreased to ^ 50% for A/[/ = 1.70 (orange circles). B The fit of the data for 
various tilts to a damped sine yields a resonance in frequency at A ?^ [/, where the amplitude re- 
verses sign (see right inset). The solid lines show the theoretically expected frequencies for the 
given parameters where the shading reflects the amplitude of the corresponding Fourier com- 
ponent. The left inset shows the three lowest eigenenergies of the Bose-Hubbard Hamiltonian 
versus the potential bias. 
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Extended two-site Bose-Hubbard model. The inclusion of next-neighbor interactions to 
the simple BHM Hamiltonian defined in Eq. SI leads to the extended two-site Bose-Hubbard 
Hamiltonian (SI) 

^EBHM ^ ^BHM _ A J ^ (n,L + h^r) {a^a^'K + al^a,,^) 

where the parameters A J and [/lr are defined in the main text. Here, we focus on a system 
of two coupled wells occupied by exactly two atoms in the two different spin-states | ]) and 
II). The relation of A J and C/lr to the tunneling matrix element J is plotted in Fig. SI. It is 
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apparent that those terms in Eq. S 1 proportional to A J act in the same way on an arbitrary Fock 
state in the system as the tunnehng operator does and therefore modify J to now ^ J + AJ, 
while the terms proportional to [/lr lead to an energy shift of the states with exactly one atom 
on each site with respect to the states with double occupancy in a single well. The states 1^) and 
|— ) stay eigenstates of the system and the states \t) and |+) are now coupled via the matrix 



which correspond to the Fourier components for the evolution of the initial state ||, |) under the 
Hamiltonian Eq. SI. In the regime of strong interactions, the frequency 002 is assigned to the 
superexchange process and by means of perturbation theory up to second order one finds the 
effective coupling parameter to be J^^ = 2J'^/C/ — t/LR. The direct spin-exchange term C/lr 
favors an antiferromagnetic ground state for repulsively interacting bosons. However, the inset 
in Fig. SI shows, that the case C/lr > 2 J'^/C/ is never reached and therefore the presence of the 
direct nearest neighbor interactions does not lead to a violation of the Lieb-Mattis theorem (S2). 
With V^ong = 10 Er and Krans = 4:1 Er, thc dependency of the two contributions on the short- 
lattice depth is well approximated by two exponential laws of the type Sg^d x exp(— |l/short /Er), 
where = 15.24£^^ for the superexchange term and = 0.3QEr for the direct exchange term. 

Loading sequence. Since the short lattice is created by a Gaussian laser beam at a wavelength 
of 765 nm, it is blue-detuned with respect to the rubidium Dl and D2 transition and thus creates 
an anti-confining potential in transversal direction. The simultaneous ramp up of the short lattice 
25 Er) together with the red-detuned transversal lattices (840 nm laser wavelength, ^ 41 Er) 
leads to a spreading of the atom distribution before freezing out tunneling and therefore to a 
more homogeneous filling than would be obtained by a fully red-detuned 3D optical lattice 
with the same beam waist and atom number. Subsequently, the long lattice is ramped up to 
^ 10 Er to obtain isolated double wells. The lattice depths during the loading and the shape of 




(S2) 



with the eigenvalues 




(S3) 
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the ramps have been optimized to give mean occupation of each double well sites close to two. 
By means of an interferometric measurement similar to the one described in (S3), we verify that 
about 70-80% of the atoms are loaded to doubly occupied sites which is compatible with the 
conversion efficiency we observe in the spin-changing collisions. 

Detection and filtering sequence. For each double well, we transfer the population of the 
left well into the second excited band of the underlying long lattice well by applying a tilt and 
ramping down the barrier afterwards (S3,S4). A subsequent adiabatic band mapping (55,5(5) 
results in spatial separation of the atoms initially localized on either side of the double well after 
TOR A Stem-Gerlach filter at the beginning of TOP allows for a spin-state selective detection. 
In addition, we transfer those atoms that have not undergone the spin-changing collisions and 
therefore stayed in |F = l,mF = 0) with a microwave 7r-pulse to |F = 2,mF = 0). They 
are removed by a resonant laser pulse during time-of-flight expansion to obtain the signal only 
from atoms in doubly occupied sites. 

Since the wave functions of both atoms are brought to overlap during the detection se- 
quence, a direct exchange term emerges which leads to spin-oscillations between the atoms 
initially located on separate sites (57). To minimize this effect, we carry out the band mapping 
immediately (~ 70 /is) after the transfer. Since we cannot suppress the direct spin-exchange os- 
cillations completely in our sequence, we believe these to limit our detection efficiency, while 
the fidelity of the preparation is probably higher and therefore closer to the expected 99%. 

Validity of the extended two-site Bose-Hubbard model. We investigate the validity of the 
EBHM by means of the frequency difference uji — uj2 in Fig. S2. While for low barrier depth, 
the EBHM models the data much better than the simple BHM, we still find deviations from its 
predictions which exceed the 2% uncertainty in our lattice depth. One reason for this can be the 
fact that the parameters for the BHM as well as for the EBHM are obtained by a single particle 
band structure calculation which does not account for the atom-atom interaction in the pairs. 
For a repulsive interaction and atom numbers larger than 1, a broadening of the single particle 
wave function has to be taken into account, which leads to lower onsite and higher intersite 
interaction energies (58,59). The comparison of the measurements to the solution of the two- 
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band Schradinger equation which includes the interaction from the start already shows a much 
better agreement. When more bands are taken into account, the two-frequency approximation 
used to fit the data is not valid anymore a direct comparison to the fit results is therefore not 
reasonable. 
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Figure SI: Extended Bose-Hubbard parameters. The parameters [/lr (light blue curve) and 
A J (black dashed curve) are plotted versus l/short- The inset shows a comparison of the direct 
exchange term [/lr (light blue) to the superexchange term 2 (dark blue), where the latter 
is always at least one order of magnitude larger than the first. 
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Figure S2: Comparison of the different models. The difference of the fitted frequencies 
a;i^2/(27r) versus l^short are plotted as in Fig. 4A and compared to the predictions of the sim- 
ple BHM (grey region, U /h), the EBHM (light green region, [U + 3[/lr)/^) and the solution 
of the two-band Schrodinger equation (dark green region). The width of the regions reflects the 
2% uncertainty in lattice depths and the error bars on the data points denote the 95% confidence 
interval as obtained from the fits. 
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